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Abstract — This paper discusses bit-level soft decoding of triple- 
parity Reed-Solomon (RS) codes through automorphism permu- 
tation. A new method for identifying the automorphism groups 
of RS binary images is first developed. The new algorithm 
runs effectively, and can handle more RS codes and capture 
more automorphism groups than the existing ones. Utilizing 
the automorphism results, a new bit-level soft-decision decoding 
algorithm is subsequently developed for general (n, n — 3, 4) 
RS codes. Simulation on (31,28,4) RS codes demonstrates an 
impressive gain of more than 1 dB at the bit error rate of 10 5 
over the existing algorithms. 

Index Terms — Reed Solomon codes, automorphism groups, 
permutation decoding, binary images, soft decoding. 



I. Introduction 

Reed-Solomon (RS) codes, with their renowned Berlekamp- 
Massey and Forney algorithms, boast robust error correcting 
capability against bit-flip errors such as large amplitude fluc- 
tuation errors and burst errors. However, to fully harness their 
power on additive white Gaussian noise (AWGN) channels 
requires effective soft decoding, which has been a research 
focus in recent years (see, for example, |]2], fl8]). 

Some algorithms, such as sorting- and scheduling- based 
message-passing algorithms, are designed for general RS 
codes 12, |H). Not restricting the underlying code specifica- 
tions, these algorithms achieve generality at the cost of a rather 
high level of computational complexity. In comparison, the 
algorithms that target specific classes of RS codes can explore 
specific code structures to effectively reduce the complexity 
and/or improve the performance. Among the latter type, a 
notable example is the bit-level automorphism-based decoder 
developed in [3|, which is designed for double-parity RS 
codes. Working especially well for short-length (double-parity) 
RS codes, the algorithm in is shown to be capable of a 
performance within 0.3dB to the maximum likelihood decoder 
(MLD) at a bit error rate (BER) of 10~ 5 . Despite its efficiency, 
however, this algorithm cannot handle RS codes of more than 
two parity symbols, in part due to its incapability of obtaining 
their automorphism groups. 

This paper studies bit-level soft decoding of triple-parity RS 
codes. Motivated by the algorithm in J3], we first investigate 
the automorphism groups of the binary images of a general 
(n, n — 3, 4) RS code, and propose a heuristic search permuta- 
tion algorithm to identify the automorphism groups. Through 
concrete examples, we show that our algorithm works effec- 
tively, can handle codes that previous algorithms cannot, and 
find automorphism groups that previous algorithms cannot. 



Exploiting these automorphism groups and integrating them 
with the general-purpose soft decoder (rather than the hard de- 
coder in O), we further develop a new bit-level soft decoding 
algorithm, termed permutation sum-product algorithm (PSPA), 
for triple-parity RS codes. Simulation demonstrates that our 
algorithm can noticeably outperform the existing algorithms, 
such as conventional hard-decision decoding (HDD) [1] and 
the soft-decision SPA |6|. The proposed new algorithm is most 
efficient for short-length RS codes, and can be extended to soft 
decode concatenated or compound codes that use short-length 
triple-parity RS codes as the component code. 

The rest of the paper is organized as follows. Section II 
introduces the binary image of RS codes. Section III discusses 
automorphism groups and details our searching algorithms for 
(n, n — 3,4) RS binary images. Section IV presents a new 
soft decoding method that effectively combines automorphism 
permutations and soft-decision sum product algorithm (SPA). 
Section V demonstrates simulation results, and Section VI 
concludes the paper. 



II. Binary Images of RS Codes 

] be a basis of ¥2^ over F2. Let the 
1. The binary image of a double- 
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[cj,o, Cj,!, Ci,2, • • • , Ci,„_i] for all j 6 Z„ and 1 < i < to. 

The parity check matrix of the (n, n— 2, 3) RS binary image 
with zeros {1, a} can be presented by a 2to x to polynomial 
in the ring F 2 [x]/(x n - 1) 0: 
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where e = a -1 , 9 £ (x) is known as the idempotent [4|, and 
9\ (x) = 1 + x + x 2 + ■ ■ ■ + x n ~ l . Examples of vectors u = 
, u m ) S Z™ are listed in Table |T] 



(ui,u 2 , 



TABLE I 
U VECTORS COMPUTED FOR 7 : 



[71:72,- ■■ ,7m]- 





Vector u 


Primitive element a 


F 23 


[2,1,0]' 


a J = a + 1 


F 2 4 


[2,1,0, 14] 1 


a 4 = a + 1 




[30,29,28,27,26]' 


a 5 = + 1 


F 26 


[4,3,2,1,0,62]' 


a b = a + 1 



III. Automorphisms of Triple-Parity RS Codes 

Previous work has reported all the automorphism groups for 
the (7,4,4) RS code over F 2 3 and part of the automorphism 
groups for the (15,12,4) RS code over F 2 4 Q. Here, we 
focus on general (n, n — 3, 4) RS codes with zeros {l, a, a 2 } 
over F 2 ™ for all rn > 4. Since, a and a 2 belong to the same 
cyclotomic coset, the parity check matrix of (n, n — 3,4) RS 
codes is represented by the following 3m x m polynomial 
matrix in the ring F2[x]/(a;™ — 1). 
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where 9(x) is short for 9 e {x). The vectors u*- 1 ' = 
[u<P,4\:.,*®\ and = [4 2) ,4 2 \... , U ^] G F™ 
can be computed using the method introduced in [13 1, and the 
results are listed in Table [TT] 



TABLE II 

Common vectors u^ 1 ' and u< 2 ) for (n, n — 3,4) RS codes. 





Vector 


Vector u^ 2 ' 




[2,1,0]' 


[2,5, 1]' 


F 2 4 


[2,1,0, 14]' 


[2,9, 1,8]' 


F 25 


[30,29,28,27,26]' 


[30,14,29, 13,28] 1 


F 2 6 


[4,3,2,1,0,62] 1 


[4, 35, 3, 34, 2, 33] 1 


F 2 7 


[6,5,4,3,2,1,0]' 


[6,69,5,68,4,67, 3] ' 


F 2 8 


[4,3,2,1,0,254,253,252]' 


[4,131,3, 130,2, 129,1,128] 1 


F 29 


[510, 509, 508, 507, 506, 


[510,254,509, 253,508, 


505,504, 503, 502] T 


252,507, 251, 506] T 


F 2 io 


[6,5,4,3,2, 1,0, 
1022, 1021,1020] T 


[6,517,5,516,4,515, 
3,514,2,513] T 



Define the symmetric group on the set {1, 2, • • • , m} as Q m . 
A permutation a e f2 m is represented using the cycle notation 
(e.g., the permutation a — (1, 2, 4) means er(l) = 2, er(2) = 4, 
and cr(4) = 1). Let id denote the identity permutation on m 
indices {1, 2, ■ • • , m}. Let a±, a 2 , ■ ■ ■ , a m g Z„ and I G Z m . 



The parameter (a,(a 1 ,a 2 , - ■ ■ ,a m ),l) completely describes 
the permutation p, such that [i, j] M> [<r(i),j2 l + a,j\. 

We present the basic idea of identifying the automorphism 
through the following example. 

Example 1: [(15, 12, 4) RS code] Let a be a fixed primitive 
element satisfying a 4 = a + 1. Consider the (15,12,4) 
RS binary image with zeros {l,a,a 2 }, imaged under the 
canonical basis. 

Consider a specific permutation p that is described by the 
parameters a — (1,2)(3,4), (01,02,03,04) and I. It can 
be verified that the permutations (cr, (a%, 0,2, ••• , a m ), I) and 
(a, (cii, <Z2j • ■ • , 0"w) + a , belong to the same automorphism 
group. Therefore, it is sufficient to fix parameter ai = 0, then 
p maps the indices 



[l,j]^[2,j2 l ], 
[3,j]^[i,j2 l +a 3 ], 



[2,j]^[l,j2 z 
[4,j]^[3,j2' 
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x)x t " 21 
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9{x)x b2i 


x)x b31 


9{x)x b32 





9(x)x b3i 


x)x bil 


9{x)x b ^ 


9{x)x bi3 






The parity check matrix of the (15,12,4) RS binary image is 
invariant under the permutation specified in (2), if and only if 
the following generator matrix in the ring F2[x]/(a; 15 — 1) is 
invariant under permutation specified in (0): 



(5) 



where (612,613,614) = (7,12,1), (631,623,634) = (8,14,4), 
(631,632,634) = (14,0,6), and (641,642,643) = (4,6,7). 

Note that permuting 9{x)x r by j n> 2 l j is equivalent to 
(9(x)x r ) = 9{x)x r , since 9{x) is an idempotent and 
9{x) 2 = 9{x). Consider [0, 9(x)x bl2 , 9(x)x bl3 , 9{x)x bli ] from 
the F2-rowspace of ((5). Applying the permutation (@), this 
row permutes to 
[0, 9(x)x bl2 , 9(x)x bi3 , 9(x)x bli ] 

[9{x)x 2lb ^+ a2 ,0,9{x)x 2,b ^+ a \9{x)x 2,b ^+ a3 ]. 
Permutation p belong to the code automor- 
phism group if the following two vector 
[9(x)x 2lb ^+ a2 ,0,9(x)x 2,b ^ +ai ,9(x)x 2,bl3+a3 ] and 
[9(x)x b21 ,Q,9(x)x b23 ,9(x)x b24 ] are generated from the 
same basis vector. This leads to the following equation 




mod (2'6i2 + a 2 - 6 CT (i) ;CT ( 2 ), n) 
mod (2'6i 3 + o 3 - 6 CT (i), CT (3),n) 
mod (2'6 M + o 4 - 6 CT (i), CT (4),") 
A x = A 2 = A 3 



(6) 



Let Ai = A 2 = A, then a 4 = mod (A + 6 



<7(l),<r(i) 



■2 l b u ,n) 



It can be seen that when Z = 0, we obtain (02,03,04) = 
(12,3,9). To verify p belong to the code automorphism group, 
we consider [9(x)x 8+r , 0, 9(x)x m+r , 9{x)x i+r ] from the F 2 - 
rowspace of (0. The permutation (@) gives rise to 
[9(x)x 8+r 7 0, 9(x)x li+r , 9(x)x i+r ] -> 

[0, 9(x)x 8+r , 9(x)x 13+r , 9{x)x 2+r ] 
which again lies in the F2-rowspace of (O. Thus, p belongs to 
the automorphism group of the (15, 12,4) RS binary images 
generated by the matrix in (|3). □ 
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To find the other permutations of a general (n, n — 3, 4) RS 
code, we propose an effective search permutation algorithm, 
which is the generalization of the above example. Note that 
there exists a set of constants {bij : 1 < i,j < m, i ^ j} in 
Z n , whereby the rows of the following in x m matrix M: 



M := 





6(x)x b2 
6{x)x b ™ 



6{x)x bl 




0(x)x b 



6(x)x bl 
6(x)x b2 





(7) 

lie in the F2-rowspace of (0. Specifically, © in the previous 
subsection is but a particular case of (Q with m = 4 and was 
reported in J3]. The proposed search algorithm is described 
below and has a complexity 0(rn\n 2 m). 

Heuristic permutation search algorithm 

Input: Polynomial matrix M in (0. 
Output: The automorphism group P. 
1. forall a e il m , 0,2,0,3 G Z„ and I 6 Z m 

Ai = mod (2 z 6i2 + a 2 - 6 CT (i) )CT (2), n) 
A 2 = mod (2 z 6i3 + 03 - ^(1)^(3), n) 
If Ax = A 2 , then a, = mod (A x + 6 CT (i),<r(i) - 2'6 li; n) 
where 3 < i < m; otherwise continue. 
Construct a permutation p = (<j, (0, 02, • • • , a TO ), Z). 
If M is invariant under p, then store p in P. 



Compute 



end 



Example 1 cont'd: [(15, 12,4) RS code] Table HIT1 presents 
all the permutations of the (15, 12,4) RS binary images with 
zeros {l,a, a 2 }, obtained by the proposed algorithm. We 
report a total of 120 unique code automorphisms, 15 of 
which map [i, j] 1-+ [i,j + a] for all a E Z15. In general, 
it appears that the permutation p = (a, (ai,a 2 , • • • ,a m ),0), 
where a — (l,n)(2, n — 1) ■ • • and (ai, <Z2>'"' , a m) = 



(0, (T 



1) 



3,6,(2- 



1 



6, 



belongs to the automorphism group of the (n, n — 3,4) RS 
binary images for all m and a 6 Z n . □ 



TABLE III 

Automorphisms of (15, 12, 4) RS code with zeros {1, a, a 2 }. 



a 


(01,02,03,04) 


I 


id 


(0,0, 0,0) + a 





(1,2)(3,4) 


(0,12,3,9) +a 





(1.3)(2,41 


(0,3, 6, 3) + a 





(MX2,3) 


(0,9,3,12) +a 





(2,4) 


(0,3, 9, 3) + a 


2 


(4,3,2,1) 


(0,0, 12, 12) + a 


2 


(1,3) 


(0,6, 0,6) + a 


2 


(1,2,3,4) 


(0,12,12,0) + a 


2 



Example 2: [(31,28,4) RS code] Consider the (31,28,4) 
RS binary image with zeros {l,a,a 2 }, where a is a fixed 
primitive element satisfying a 5 — a 2 + 1, imaged under the 
canonical basis. The parity check matrix of this binary image 
is invariant under permutation p, if and only if the following 
polynomial matrix in the ring F2[x]/(x 31 — 1) is invariant 



under permutation p: 

8(x)x 6 
0(x)x 22 
0(x)x u 
9{x)x 
9{x)x 28 



6(x)x 22 
9{x)x 12 
e{x)x 2 



6(x)x 15 
9{x)x 11 
9(x)x 5 


6{x)x 21 



6(x)x 2A 
0(x)x 14 
9(x)x 11 
6(x)x 3 




(8) 



9{x)x 12 
0(x)x 5 


8{x)x 20 
6{x)x 13 

Applying the proposed search algorithm leads to Table |IV] 
We identify a total of 124 unique code automorphisms, 31 of 
which map +» + a] for all a G Z31. □ 

TABLE IV 

Automorphisms of (31, 28, 4) RS code with zeros {l, a, a 2 } . 



a 


(01,02,03,04,05) 


I 


id 


(0, 0,0,0, 0) + a 





(1,2)(4, 5) 


(0,15,23,29,17) +a 





(M)(2,5) 


(0,29,9, 18,20) + a 





(1,5)(2,4) 


(0, 17,3, 20,6) + a 






Example 3: [General (n, n — 3, 4) RS codes] The number 
of all the automorphism groups found via the proposed algo- 
rithm is listed in Table [V] for several (n, n — 3, 4) RS codes. 
We also present the code automorphism orders computed using 
MAGMA software. Compared to the existing algorithms, 
our algorithm finds significantly more automorphisms than 
the method in J3]; It produces results consistent with those 
obtained from the MAGMA computation in [5Q, but can 
handle fields of larger orders (e.g. F 2 s) where MAGMA would 
face difficulty. □ 

TABLE V 

Automorphism subgroup orders of (n, n — 3, 4) RS codes over 

FINITE FIELD F 2 m OBTAINED BY THE PROPOSED SEARCH ALGORITHM. 
(N.A. MEANS NOT AVAILABLE) 





Proposed 


MAGMA |5| 


Lim et at | 3 1 


F 2 4 


120 


120 


18 


F 2 5 


124 


124 


N.A 


F 26 


126 


126 


N.A 


F 2 7 


254 


254 


N.A 


F 2 8 


510 


N.A 


N.A 



IV. Soft Decoding of (n, n - 3, 4) RS Codes 

The availability of a list of the permutation groups allows 
us to develop more effective soft decoders for triple-parity RS 
codes. The proposed new permutation sum-product algorithm, 
illustrated in Fig. Q] takes advantage of the rich permutation 
groups obtained in Section [III] and the soft decision output of 
the classical (binary) sum-product algorithm. 



Received 


®* 


Sum-Product 


® 




codeword y 




decoder | 







Random 
permute 

v, =p(y) 



© 



y+r>'(y„) 



© 



Corrected 1 © 


Hard decision 


codeword y J 


decoder 



Fig. 1. Permutation sum-product algorithm for (n,n — 3,4) RS binary 
images. 

'Note that two groups of the same order do not imply the same. 
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Permutation sum-product algorithm (FSFA) 
Input: Observation y 

Output: corrected codeword. 

1. Represent y as a binary image. Decode y using the 
conventional sum-product algorithm (SPA). 

2. If the SPA produces a valid codeword c, then stop; otherwise 
store the log-likelihood y 

3. Generate a random permutation y: y p = p(y) (e.g. using 
the proposed search algorithm). 

4. Repeat step 1 until all the permutations are exhausted or a 
maximum number of trials is reached. 

5. Compute the sum: y sum = y + £ p _1 (y p ). 

6. Perform hard decision on y sum and outputs the result. 

Remark: It should be noted that our decoding algorithm 
differs from that proposed in [0. The algorithm in |]7] is 
restricted to binary cyclic codes over binary erasure channels. 
Our problem here involves the general permutations of bit- 
level RS codes over AWGN channels, and the cyclic codes in 
may be viewed as a special case of ours. Our decoder also 
differs from that in [3|. To achieve its desirable performance, 
the decoder in [3 1 must be able to find a "perfect" permutation 
that gathers all of the bit errors in one symbol, and then decode 
this permutation using some hard-decision decoder. Not only 
is such a perfect permutation very difficult to identify, but the 
decoder performance would also degrade very quickly should 
a perfect permutation become unavailable. In comparison, our 
PSPA method obviates the burden of identifying a perfect per- 
mutation, simply takes in random permutations, and exploits 
the power of the soft-decision SPA to effectively correct bit 
errors. 

V. Simulation results 

Simulations are conducted to verify the effectiveness of 
the proposed decoding algorithm. The triple-parity (31,28,4) 
and (63,60,4) RS codes over F 2 s and F 2 e are considered. 
We compare our proposed soft-decision PSPA algorithm with 
three other systems: the uncoded system, the RS code with 
the conventional hard-decision decoding (HDD) implementing 
the Berlekamp-Massey algorithm (TJ, and the RS code with 
the soft-decision sum-product algorithm (SPA) [6|. Since the 
decoder in [3 | only work for double-parity RS codes, it is not 
included in the comparison. 

In the proposed PSPA, the maximum number of decoding 
iteration is set to 30 for the embedded sum-product decoder, 
and the number of random permutations (trials) is set to 10 
(i.e. at the most 10 automorphism permutations are chosen 
randomly from the 4 classes in Table IV). As shown in Fig. 
|2] and Fig. [3] our PSPA noticeably outperforms the hard- 
decision decoder and the soft-decision SPA by more than 1 
dB at a BER of 10~ 5 . This impressive performance gain is 
attainable because automorphism permutations lead to new and 
different decoding matrices, opening the possibility to break up 
some critical short cycles or trapping sets that were previously 
preventing the SPA algorithm from converging. 

VI. Conclusion 

We have investigated the bit-level soft decoding of triple- 
parity (n, n — 3, 4) RS codes. We first proposed an effective 




Signal to Noise Ratio (SNR) 

Fig. 2. BER performance for (31, 28, 4) RS code over AWGN channels. 
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Fig. 3. BER performance for (63, 60, 4) RS code over AWGN channels. 

search algorithm to identify automorphism groups of a general 
(n,n — 3,4) RS code with zeros {l,a, a 2 }, and next devel- 
oped a soft-decoding algorithm that integrates automorphism 
permutation and the conventional soft-decision sum product al- 
gorithm. Simulations confirm the effectiveness of the proposed 
permutation sum-product algorithm, revealing an impressive 
performance of more than 1 dB at the BER of 10~ 5 than the 
existing algorithms. Our algorithm is particularly suitable for 
short-length triple-parity RS codes, or concatenated codes in 
which they are a component code. 
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